On the transverse mode of an atom laser 
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The transverse mode of an atom laser beam that is outcoupled from a Bose-Einstein condensate is 
investigated and is found to be strongly determined by the mean-field interaction of the laser beam 
with the condensate. Since for repulsive interactions the geometry of the coupling scheme resembles 
an interferometer in momentum space, the beam is found show filamentation. Observation of this 
effect would prove the transverse coherence of an atom laser beam. 
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Research on atom lasers is an active and fascinating 
area in atomic physics Q [21 • Several laboratories around 
the world are now using continuous output couplers to 
produce atom laser beams from Bose-Einstein conden- 
sates. It is therefore important to characterize the qual- 
ities of these beams. Recently, their temporal coherence 
was verified Q and their transverse divergence was mea- 
sured 01 . 

Although atom lasers and optical lasers show strong 
similarities, the possibility for atoms to scatter off each 
other leads to various effects absent in optical lasers. In 
the following we will show that a strong, inhomogcneous 
repulsive potential, as is represented by the remaining 
Bose-Einstein condensate, can be a source of instabilities 
for the beam and, in particular, can lead to its transverse 
filamentation. Since most experimental setups involve 
asymmetric traps with two stiff directions, our results 
are immediately applicable to these experiments. 

In magnetic traps output couplers for atom laser beams 
are realized by coupling a fraction of the Bose-Einstein 
condensate into a magnetically untrapped state. This 
process happens inside the trapped sample along a sur- 
face where the resonance condition for output coupling 
is fulfilled, subjecting the output coupled atoms to the 
repulsive mean field potential of the condensate. Since 
gravity displaces the symmetry axis of the Bose-Einstein 
condensate with respect to the symmetry axis of the mag- 
netic trapping field, the repulsive potential is however 
not homogeneous over the resonance surface. In a classi- 
cal picture this situation corresponds to particles rolling 
off a potential from different heights, leading to a non- 
negligible momentum spread or dispersion in the trans- 
verse directions Moreover, the Snite interaction time 
of the falling beam with the remaining condensate leads 
to a non-monotonic increase of the transverse atomic po- 
sition x{t) as a function of the points of resonant output 
coupling Xi. In the quantum dynamics, this leads to the 
interference of atoms with different transverse momenta 
within the beam. 

In the following we will investigate this interference 
process. First we consider an idealized, one-dimensional 
model in the direction perpendicular to gravity, and 



we examine its classical, semiclassical and quantum be- 
haviour. After we have identified and described the rele- 
vant processes, we will include the effects of gravity and 
propagation in more than one dimension. 

Atom lasers with radiofrequency output couplers usu- 
ally couple different Zeeman substates of the trapped 
atoms. For a Bose-Einstein condensate in a |F = 
l,mF = —1) state, two output states are possible. Ei- 
ther mp ~ 01 mp = 1, the first of which has a van- 
ishing magnetic moment and the latter experiences a re- 
pulsive force by the magnetic trap. However, since the 
output coupling rate from the mp = — I state into the 
mp = state is usually chosen to be small, subsequent 
transitions into the mp = 1 state can be neglected. We 
therefore restrict our considerations to a two-level sys- 
tem, where the important coupling parameters are the 
Rabi frequency H. and the detuning A of the rf-field. 
For weak coupling the resonance condition is determined 
by the spatial dependence of A within the Bose-Einstein 
condensate. After a transformation into a co-rotating 
frame tpmp (t) e"*'"'"'^'^* Vm^ (i) followed by the stan- 
dard rotating wave-approximation, the equations for the 
condensate wave function ipc and the atom laser beam 
wave function tpi, are given by a set of coupled Gross- 
Pitaevskii-equations [g 

+^/(|V^P + IV'Jf )V'. + ?if^V'j- , 

with i,j = c,h and U = ^Trh'^as/m. We have as- 
sumed that all triplet scattering lengths have the same 
value ttij = as and we will choose Os to be pos- 
itive. The external potentials are given by T4(r) = 
Y {yj\{x'^ + z^) + i^yV'^) + mgz and Vfc(r) — mgz, where 
g is the gravitational constant, luj_ and ujy are the trap- 
ping frequencies of the cylindrically symmetric magnetic 
field and m denotes the mass of the atoms. 

It has been found, using a separation ansatz for the 
spatial modes, that the atom laser beam in the direc- 
tion of gravity can be almost perfectly described by an 
Airy-Function H H [13 . We therefore first consider the 
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behavior in the transversal x-direction as independent of 
the other directions. 

A naive, strictly one-dimensional treatment of the 
atom laser in the horizontal x-direction would result in 
exactly two resonance points x = ± ^ 2hujrf / muj'^ 
This, however, is not a good approximation to the three 
dimensional situation of the experiment. Since gravity 
leads to a displacement of the condensate from the cen- 



ter of the magnetic field by an amount 



the resonance shell crosses the condensate with very low 
curvature and is therefore better approximated by a 
plane. This means that in the horizontal directions out- 
put coupling happens along the full Thomas-Fermi dis- 
tribution and the initial beam wave function is a scaled 
down copy of the condensate wave function 

While the atoms fall under gravity, the mean field po- 
tential they experience from the condensate changes. To 
account for this in the one-dimensional approximation 
one has to diminish the mean field potential, Uc{x,t) — 
U\'4>c{x,zt)\'^ during the evolution according to the free 
fall of the atoms, zt = Zg + 

However, the principal physics of the horizontal mode 
and its instability is most clearly demonstrated by first 
considering the case with g = 0, i.e. taking Uc{x,t) = 
Uc{x,Zg) to be constant in time. The effective one- 
dimensional Gross-Pitaevskii equation for the beam is 
then given by 



ih, 



dt 



2m dx^ 



+ Uc{x) + U\A\ 



(2) 



In the Thomas-Fermi approximation the mean field po- 
tential is given by a truncated inverted harmonic oscilla- 
tor 



Uc{x) = 



11 [l- 




for < xtf, 
for > xtf- 



(3) 



Here /i is the chemical potential and xtf = \/2[ilmu)\ 
the Thomas-Fermi radius of the condensate. 

Let us first analyse eq. lj2Jl classically and neglect the 
nonlinear term t7|^f,p, because it is normally three or- 
ders of magnitude smaller than Uc{x). Neglecting for 
a moment the truncation of the potential at the con- 
densate border, i.e., assuming the potential Uc{x) = 
— x^ /x'^p),\/x, the classical equation of motion 



dUcjx) 
dx 



(4) 



can be exactly integrated by x{t) = Xi coah{uj±t), with 
Xi the initial position of the atoms at the time of the 
outcoupling. Since the cosh is an exponentially increasing 
function of the time, the resulting evolution is a spreading 
of the initial distribution. This means that atoms with 
larger Xi will also have larger x at later times t, (see 
dashed fine in Fig.^). 

The density distribution using the potential in eq. Q is 
shown in Fig. |21 (dashed lines). The atoms show a tempo- 
rary localization at the Thomas-Fermi radius, since the 
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FIG. 1: (a) Position of atoms starting at Xi{t = 0) < xtf 
and rolling off an inverted harmonic oscillator potential for 
a time of f = 2 ms. The full line shows the results for a 
potential that is truncated at xtf ~ 4^m (cf. eq. (O) and the 
dashed line shows the results for an untruncated potential, (b) 
Time needed for the atoms starting at a = Xi/xTF to reach 
the condensate boundary at xtf (cf. eq. The chemical 

potential in both calculations \s = 1700Hz. 



arrival time distribution for the atoms at the Thomas- 
Fermi edge 



tTFict) — — hi I — 

uj \ \ a 



(5) 



with a = Xi/xTF, shows a plateau (see Fig. ^). Once 
the atoms have passed the Thomas-Fermi radius all their 
initial potential energy, E = ^(1 — a^), has been trans- 
formed into kinetic energy and atoms starting closer to 
the condensate center {x ~ 0) therefore end up with 
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FIG. 2: Transverse density distributions of the beam obtained 
from a classical (dashed) and a quantum-mechanical calcula- 
tion for two different evolution times. The figures in the lower 
row show the absolute square of the Fourier transform of the 
quantum mechanical wave function. A strong filamentation 
is clearly visible in the quantum results. 
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higher final velocity. In the asymptotic limit, atoms origi- 
nating from the center will have overtaken all other atoms 
and the density distribution that originally had a nega- 
tive slope will have a positive slope. This can be seen 
from the dashed curves in Fig. |21 

The full quantum mechanical behaviour can be found 
by solving Eqs.JQ) numerically, for which we use a stan- 
dard split-operator/FFT technique. As can be seen from 
Fig. m (solid lines), although the general feature of local- 
ization is preserved, the density distribution is modulated 
by an interference pattern. The reason is that for finite 
times, t > txF, atoms with different initial positions ar- 
rive at X at the same time t, because the equation 
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XTF 



(l + ujy/l - {t - tTF(a))) 



(6) 



has more than one solution for a, and quantum mechan- 
ically, atoms coherently outcoupled at different a's in- 
terfere. Observation of these fringes would prove the 
transversal coherence of the atom laser beam. 

A mathematically and conceptually very elegant 
method to calculate the interference pattern is by use 
of path integrals p!3| • The unitary wave function prop- 
agator is obtained by adding phase factor contributions 
over all paths x{t) along which an argument Xi in the 
initial wave function can evolve into the argument x j of 
the final state 



(7) 



The phase factor is the classical action along the path 

1 



s[x{t)] = y dt 



—rax 
2 



V[x{t)] 



(8) 



Path integrals are usually not exactly solvable, but they 
can be approximated well when for example few classical 
paths Xci{t) are dominant. We can then restrict the inte- 
gral to small deviations around them S[x\ = S[xci + 5x]. 
This will provide a qualitative understanding of the ob- 
served fringes as well as a very good account of the quan- 
titative results of the quantum calculations. 

From Fig.^one can immediately make the observation 
that for x > xtf two classical paths contribute. The 
atoms emerge from two different initial positions, and 
they will therefore arrive at x with two different momenta 
fiki and hk2. This suggests interference fringes of width 
(|fci — fc2|)~^, i.e., larger fringes for large values of x (see 
inset of Fig|3Jl. 

For potentials whose second spatial derivative is con- 
stant, the propagator can be written as 



dx 



k{xY 

. (9) 

with ki — k{xi), kf — k{xf) and the classical action is 
given by 



^[xcz] = —E {tf — ti) + h dx k{x) 



(10) 
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FIG. 3: Transverse density distribution for t=5ms of the 
beam from the quantum-mechanical calculation (solid line) 
and from the analytical expression based on Feynman path in- 
tegrals as explained in the text (dashed line) . The inset shows 
the distance 5x between successive minima imin of |V'P, start- 
ing at the Thomas-Fermi radius. The units of x is meters. 
The disagreement in amplitude stems from different initial 
conditions due to the different numerical methods used. 



Inserting the inverted and truncated harmonic oscillator 
potential of eq. V{x) — Udx), the action can be 
calculated to give 



S[Xc 



9 2xf ~ Xtf I — 
-a^i^itTF H V 1 



XTF 



(11) 



for Xj > Xtf and the bracketed expression in eq. lO can 
be calculated to be 



2TTikfki 



dx 
k{xy- 



1 



a 



(12) 



where oo = yjh/muj± . In Fig. |3 the evolution of the 
wave function according to the semiclassical approxima- 
tion to eq. (O is shown (dashed line). Comparison with 
the full quantum mechanical evolution of eq. Q shows 
almost perfect agreement. 

To justify and understand the semiclassical approxi- 
mation, let us note that restriction to classical paths and 
the assumption of quadratic potentials is formally equiv- 
alent to a WKB approximation and the condition of 
validity can be written as 



1 



dk 
dx 



1 



< 1 



(13) 



where we have used the momentum of the atoms given 
by hk{x) — ^2m[E — Vix)). This condition is well ful- 
filled as long as V(x) is an inverted harmonic oscillator 
for all atoms. One may, however, question the validity 
close to the Thomas-Fermi radius, where the slope of 
the potential changes strongly within the healing length 
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FIG. 4: Two dimensional simulations of the atom laser, (a) 
The beam after an evolution time of 2.5 ms. (b) A cut through 
the density distribution ai z = — 15/im, corresponding to 
an evolution time of 1.3 ms. (c) The corresponding Fourier 
transform. Comparison with Fig. |21 shows good agreement. 



of the system, i.e., one might expect the appearance of 
quantum reflection effects at this point. To estimate 
the reflectivity of the potential step, we note that the 
de-Broglie-wavelength of the atoms when reaching the 
Thomas-Fermi radius is larger than the healing length f 



AdB = 



> 



TF 



27rao 

XTF 
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XTF 



(14) 



One can therefore approximate the edge of the conden- 
sate by an effective step. Since the exact choice of the 
position of the effective potential step is not crucial, one 
can from eq. H14|l estimate its effective height to be as 
large as 10% of the central mean field potential. This 



leads to a reflection coefficient 

R = 



1 - 20Vl - a2 



< 10- 



(15) 



justifying very well the scmiclassical treatment. 

Let us finally consider a two dimensional situation in- 
cluding gravity. The beam atoms are subjected to the 
mean-field potential for a time of the order of 0.5 — 1ms 
assuming a typical condensate radius of 4^m. Most of 
the atoms therefore do not experience a complete hor- 
izontal roll-off from the mean-field potential (compare 
Fig.^), however qualitatively the above picture remains 
unchanged. We have solved the two-dimensional version 
of eq. ^ in the x- and z-plane numerically and in Fig.^ 
the atom laser beam for short evolution times is shown. 
Once the beam has left the overlap area with the Bose- 
Einstein condensate its evolution within the transverse 
direction is completely determined by a free evolution 
and the far-field result can be simply calculated by the 
Fourier transform. A cut through the density distribu- 
tion is shown in Fig.^ and the far field of this distribu- 
tion is shown in Fig. Both pictures show good quali- 
tative agreement with the results of the one-dimensional 
analysis. 

In summary we have shown that the transverse mode 
of an atom laser is strongly determined by the interaction 
of the beam with the mean-field of the residing conden- 
sate. Due to the finite time of this interaction, the system 
resembles an interferometer in momentum space and the 
beam shows filamentation in the transverse directions. 
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